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Abstract -Using Landau’s theory of two-fluid hydrodynamics we investigate the sound modes 
propagating in a uniform weakly-interacting superfluid Bose gas for values of temperature, up 
to the critical point. In order to evaluate the relevant thermodynamic functions needed to solve 
the hydrodynamic equations, including the temperature dependence of the superfluid density, we 
use Bogoliubov theory at low temperatures and the results of a perturbative approach based on 
Beliaev diagrammatic technique at higher temperatures. Special focus is given on the hybridization 
phenomenon between first and second sound which occurs at low temperatures of the order of the 
interaction energy and we discuss explicitly the behavior of the two sound velocities near the 
hybridization point. 
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Second sound is one of the most interesting manifes¬ 
tations of superfluidity. It is associated with the motion 
of the fluid where its normal and superfluid components 
oscillate in opposite phase (for a recent discussion see 
[1]). In superfluid Helium, where it was first measured 
by Peshkov [2], second sound is usually regarded as a tem¬ 
perature (or entropy) wave and is distinct from first sound 
which is basically a density wave. A major interest of sec¬ 
ond sound is that it gives direct access to the superfluid 
density [3], a quantity of high physical interest, vanish¬ 
ing at the superfluid critical temperature. Recently the 
measurement of second sound has been used to determine 
experimentally the temperature dependence of the super- 
fluid density in a strongly interacting atomic Fermi gas [4]. 
From a theoretical point of view the description of first and 
second sound relies on Landau’s two-fluid hydrodynamic 
equations of superfluidity [5]. Below the critical temper¬ 
ature, with the appearance of a superfluid component in 
the system, the macroscopic dynamics of the system is de¬ 
scribed by coupled equations involving the normal and the 
superfluid components. The peculiar nature of these equa¬ 
tions arises on the one hand from the hydrodynamic na¬ 
ture of the superfluid motion, strictly related to the phase 
of the order parameter, and on the other hand from the 
collisional regime characterizing the normal component. 


In the linearized limit of small amplitude oscillations, the 
two-fluid hydrodynamic equations take the form 


d 2 p 

dt? 


V 2 P 


( 1 ) 


9 2 S = Ps5 2 v 2 T 
dt 2 p n 


( 2 ) 


where p s and p n ( p s + p n = p = mn ) are the densities of 
the superfluid and normal components, s is the entropy 
density per unit mass and P is the pressure of the gas. 

From Eqs. (1-2), one derives the following quartic equa¬ 
tion for the sound velocities (see for example [6]) 


c 


4 


fdP 

p s Ts 2 8 

| c2 + PsTs 2 OP 

\dp 

g Pn,Cv / 

' PnC v dp 


0 


(3) 


where c v is the specific heat per unit mass and ( dP/dp)r 
and ( dP/dp)s are the inverse isothermal and adiabatic 
compressibilities respectively. Equation (3) admits simple 
exact solutions at low T and close to the critical tem¬ 
perature. In the first case, where all the thermodynamic 
functions are fixed by the thermal excitation of phonons, 
one can easily find the Landau’s result c 2 = cf/3 where 
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ci = \J ( dP/dp)r=o is here the T = 0 first sound veloc¬ 
ity. Close to the critical temperature, where the super¬ 
fluid density vanishes in 3D systems, the second sound 
velocity instead tends to zero. In this limit one can ne¬ 
glect, in Eq. (3), the small c 4 term as well as the sec¬ 
ond term, proportional to p Sl in the parenthesis mul¬ 
tiplying c 2 . Then, using the thermodynamic relation 
c p /c v = ( dP/dp)g/(dP/dp)T , one gets the result: 
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In the same limit the first sound velocity approaches the 
usual expression 
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In systems characterized by a small value of the thermal 
expansion coefficient, such as liquid 4 He, the adiabatic and 
isothermal compressibilities are nearly the same and equa¬ 
tion (3) admits the very simple solutions c 2 = ( dP/dp ) 
and c 2 = ( p s Ts 2 )/(p n c ) holding for all temperatures be¬ 
low T c where the specific heat c can be evaluated at either 
constant volume or pressure. 

A major question concerns the behavior of the sec¬ 
ond sound velocity at intermediate temperatures in dilute 
gases. In superfluid Fermi gases with resonant interactions 
the thermal expansion coefficient is still relatively small 
and the formula (4) for the second sound velocity turns 
out to be accurate in a wide temperature interval. At 
the same time the thermal expansion is still large enough 
to show up in a detectable signal in the density fluctu¬ 
ations associated with the propagation of second sound. 
Thanks to this coupling it has been actually possible to 
measure the velocity of second sound in Fermi gases inter¬ 
acting with large values of the scattering length in highly 
elongated configurations and to extract the temperature 
dependence of the superfluid density [4,7]. 

Weakly-interacting Bose gases behave in a very differ¬ 
ent way due to their large isothermal compressibility. An 
important consequence of this different behavior is the oc¬ 
currence of a hybridization phenomenon between first and 
second sound [8-11]. The mechanism is caused by the 
tendency of the velocity of the two modes to cross at very 
low temperatures, of the order of the T = 0 value of the 
chemical potential p(T — 0) = gn where g — dirh 2 a/m is 
the bosonic coupling constant, with the consequent emer¬ 
gence of a hybridization mechanism. For temperatures 
below the hybridization point the velocity of the upper 
branch approaches the Bogoliubov value Cb = y/gn/m, 
while the lower branch approaches the Landau’s result 
C 2 = \Jgn/Srn. Above the hybridization point the roles of 
first and second sound are inverted, in the sense that the 
second sound is essentially an oscillation of the superfluid 
density which, for temperatures significantly smaller than 
the critical temperature, coincides with the density of the 
gas. 


The main purpose of the present work is to exploit 
in an explicit form the mechanism of hybridization, pro¬ 
viding an exact expression for the hybridization tem¬ 
perature in the case of weakly-interacting Bose gases 
and to discuss the behavior of second sound both be¬ 
low and above the hybridization point. For this pur¬ 
pose we will use Bogoliubov theory to calculate the rel¬ 
evant thermodynamic functions entering Eq. (3) at low 
temperatures (T <C T c ), in terms of the elementary ex¬ 
citations of the gas, whose dispersion is given by the fa¬ 
mous Bogoliubov law e(p) = y/p 2 /2m (p 2 /2m + 2 gn) and 
which are thermally excited according to the bosonic rule 
N v [e) = (e e( - p ^ kT — 1) _1 . Instead, at higher tempera¬ 
tures we will use the perturbation theory developed in 
Ref. [12] based on the Beliaev diagrammatic technique at 
finite temperature [13]. In this approach thermal effects 
in the Bogoliubov excitation spectrum are accounted for 
through a self-consistent procedure. At low temperatures 
this approach coincides with Bogoliubov theory, except 
for temperatures smaller than ( na 3 ) 1 ^ i gn , i.e. at tem¬ 
peratures much smaller than the hybridization point. At 
higher temperatures the theory of Ref. [12] turns out to 
be very accurate in dilute gases when compared with ex¬ 
act Monte-Carlo simulations. It follows that, at least for 
small values of the gas parameter na , Bogoliubov theory 
and the diagrammatic approach of Ref. [12] match ex¬ 
actly in the hybridization region of temperatures kT ~ gn 
and that the thermodynamic behavior of the gas is conse¬ 
quently under control for all ranges of temperatures, both 
below and above gn. 

The coefficients of the quartic equation (3) depend not 
only on the equilibrium thermodynamic functions, but 
also on the normal and superfluid densities. The normal 
density can be calculated using the Landau’s prescription 


Pn 


1 f dATp(g) 2 dp 
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in terms of the elementary excitations of the gas. The Lan¬ 
dau’s prescription (6) ignores interaction effects among el¬ 
ementary excitations and in a dilute Bose gas is expected 
to be very accurate except for temperatures close to the 
critical point. 

A peculiar property of the weakly-interacting Bose gas is 
that all the thermodynamic functions entering Eq. (3), as 
well as the normal density p n , can be written in a rescaled 
form as a function of the reduced temperature t = kT/gn 
and of the reduced chemical potential 77 = gn/kT// where 


kT" = 


2irh~ 


n \ 
to VC(3/2 )) 


2/3 


(7) 


is the critical temperature of the ideal Bose gas. In weakly- 
interacting Bose gases Xj? does not coincide with the ac¬ 
tual critical temperature which contains a small correction 
fixed by the value of the gas parameter na 3 : 

T c =T c °(l + 7 (na 3 ) 1 / 3 ) (8) 
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with 7 ~ 1.3 [14-16]. The reduced chemical potential 
can also be expressed in terms of the gas parameter as 
g = 2C(3/2 ) 2 / 3 (na 3 ) 1 / 3 . 

Using Bogoliubov theory, the free energy F = U — TS 
can be written as 
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(9) 

where we have defined p = p/y/mgn and, for complete¬ 
ness, we have included the Lee-Huang-Yang correction to 
the t = 0 value of the free energy (term proportional to 
(na 3 ) 1 / 2 ). 


The thermodynamic functions entering Eq. (3) are then 
easily calculated using the following thermodynamic rela¬ 
tions (p = mn) 


dP 

~d~ P 



1 dF/N 



m dT 

P 

_ 1 d(F/N + mTs ) 

'-'V 

m 

d T 

f- 

dP 

kT 

( dP 

.r 4 

+ 2~ 
rp P Cy 

dT 


( 10 ) 

( 11 ) 


( 12 ) 


where P = — (dF/dV)x is the pressure of the gas. 
The normal density ( 6 ) instead takes the form 
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2 dp 
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The idea now is to calculate the two solutions of Eq. (3) 
for a fixed value of t of the order of unity, taking the limit 
g — > 0. Physically this corresponds to considering very low 
temperatures (of the order of gn) and very small values of 
the gas parameter na 3 . For example in the case of 87 Rb 
the value of the scattering length is a = 100ao (where 
ao is the Bohr radius) and typical values of the density 
correspond to na 3 ~ 10~ b . This corresponds to g « 0.04 1 . 

Writing the solutions of Eq. (3) in terms of the T = 0 
Bogoliubov velocity c B = \J gn/m , one finds that, as 
g —> 0 , the two sound velocities only depend on the di¬ 
mensionless parameter t and are given by 


r 2 - r 2 
— C B 


c - = 4/4 


(14) 
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Fig. 1: Sound velocities (light and dark blue) computed by in¬ 
terpolating Bogoliubov theory and the diagrammatic approach 
of Ref. [12], over the whole range of temperatures 0 < t < t c . 
The gas parameter is chosen to be na 3 = 10 -fi . The critical 
point (see Eq. ( 8 )) corresponds to t c — kT c /gn = 26.7. The 
inset shows the hybridization region. 


that / is consequently a function of t, independent of g. 

In the g —> 0 limit, the two velocities shown in Fig. 1 
cross each other at the value th y b ~ 0.6. At lower temper¬ 
atures c 2 approaches, as expected, the zero-temperature 
value 4/3- 

By considering finite, although small, values of g , it is 
possible to show that, at the hybridization point, the two 
branches exhibit a gap proportional to g 3 / 4 . The mech¬ 
anism of hybridization, explicitly shown in the inset of 
Fig. 1, permits to identify an upper branch C\ (which co¬ 
incides with c+ for t < fhyb and with c_ for t > fhyb) 
called “first sound”. The lower branch C 2 (called “second 
sound”) instead coincides with c_ for t < fhyb and with 
c+ for t > fhyb- 

The validity of Eqs. (14) is limited to very low tem¬ 
peratures where the thermal depletion of the condensate 
can be ignored and Bogoliubov theory can be safely ap¬ 
plied. When the temperature starts being comparable 
with the critical temperature T c (corresponding to the 
value t c = 26.7 in the case of Fig. 1, where we have cho¬ 
sen g = 0.04), Bogoliubov theory is no longer applica¬ 
ble. In fact at such temperatures the thermal depletion 
of the condensate becomes important and the Bogoliubov 
expression for the dispersion law is inadequate. The su¬ 
perfluid density fraction, calculated according to Eq. (13) 
with the T = 0 value of the Bogoliubov dispersion law, 
vanishes at T ~ 1.2T C , well above the critical temperature 
(8), further revealing the inadequacy of the theory at high 
temperatures. As anticipated above a better approach to 
be used at temperatures of the order of the critical value is 
the diagrammatic approach of Ref. [12] whose predictions 


where /(t) = lim^o ^77 ^ ■ Using Eqs. (9-13), one ac¬ 
tually finds that s oc g 3 / 2 , c v oc g 3 ^ 2 and p n oc g 3 ^ 2 and 


1 Notice that such a value is significantly smaller than the ratio 
fi(T = 0 )/kT c evaluated in the presence of harmonic trapping with 
the same value of the gas parameter na 3 where is n is the density in 
the center of the trap [6]. 
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for the first and second sound velocities at temperatures 
higher than the hybridization temperature are reported in 
Fig. 1 and which properly interpolates with the predictions 
of Bogoliubov theory near the hybridization point. 

For temperatures significantly higher than the hy¬ 
bridization temperature an accurate and elegant expres¬ 
sion for the second sound velocity is obtained by evalu¬ 
ating all the quantities entering the quartic equation (3), 
except the isothermal compressibility and the superfluid 
density, using the ideal Bose gas model. One can actually 
show that the adiabatic compressibility, the specific heat 
at constant volume and the entropy density deviate very 
little from the ideal Bose gas predictions in a wide interval 
of temperatures above the hybridization point. This sim¬ 
plifies significantly the calculation of Equation (3). In fact 
in the ideal Bose gas model one finds the simple result 

( |L) + *Zii-r2i (15) 

\dpjj Pn C v p n C v 


for the coefficient of the c 2 term. It is now easy to derive 
the two solutions satisfying the condition C\ C 2 . To 
obtain the larger velocity ci one can neglect the last term 
in (3). Using the thermodynamic relations of the ideal 
Bose gas model and identifying the normal density with 
the thermal density ( p n = mriT = mN^/V) one 
the prediction 

2 _ 5 <75/2 ksT 

1 3g 3/2 TO 

for the first sound velocity (Lee and Yang ,1959). 
culate C 2 we must instead neglect the c 4 term in 
using result (15) one finds the useful result 


obtains 


(16) 

To cal- 
(3) and 
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(17) 


revealing that the superfluid density and the isothermal 
compressibility are the crucial parameters determining the 
value of the second sound velocity of weakly-interacting 
Bose gases for T gn. In the relevant region 10 < t < 25 
of the reduced temperature, the expression Eq. (17) for the 
second sound velocity deviates from the exact solution of 
the two-fluid hydrodynamic equation reported in Fig. 1 by 
less than 7 percents. 

In Fig. 2 we show the temperature dependence of the su¬ 
perfluid density predicted by the diagrammatic approach 
of Ref. [12] compared with the predictions of Bogoliubov 
theory, the results of quantum Monte Carlo simulations 
as well as the ideal Bose gas value p( 1 — (T/T°) 3 / 2 ). No¬ 
tice that the diagrammatic approach of Ref. [12] predicts 
an unphysical jump of the superfluid density at the tran¬ 
sition point and that critical fluctuation terms should be 
taken into account for a correct description near T c . It 
is also worth discussing the behavior of the isothermal 
compressibility (see Fig. 3), a quantity of high interest 
near phase transitions. The diagrammatic approach of 



Fig. 2: Superfluid density computed using the ideal Bose gas 
model (orange curve), the Bogoliubov theory (blue curve) and 
the diagrammatic approach of Ref. [12] (green curve). The 
black dots correspond to Monte-Carlo simulations. The pa¬ 
rameters are chosen as in Fig. 1. 


Ref. [12], which at high temperatures approaches Hartree- 
Fock theory, predicts an increase of the isothermal com¬ 
pressibility (pdP/dp\T)~ 1 with respect to its T = 0 value 
(gn)^ 1 and a divergent behavior near the critical temper¬ 
ature, which is typical feature of mean field theories. The 
inclusion of critical fluctuations will modify this behavior 
near the critical point. To our knowledge the isothermal 
compressibility near the critical point of dilute Bose gases 
has not yet been accurately calculated nor measured. This 
contrasts with resonant Fermi gases where its recent mea¬ 
surement [17] was actually employed to identify experi¬ 
mentally the value of the critical temperature. 

In Fig. 4, we finally show the ratio (8n/n)/(ST/T) be¬ 
tween the relative density and temperature variations cal¬ 
culated for the first and second sound solutions of Eq. (3) 
as a function of temperature. This quantity represents an 
important characterization of the two branches. It is in 
fact well known that sound in an ideal classical gas is an 
adiabatic oscillation characterized by the value 3/2 for the 
ratio (Sn/n)/(ST/T). For the upper solution of the hydro- 
dynamic equations (first sound), this ratio is positive over 
the full range of temperatures below the transition and its 
value increases with T, getting close to unity near the tran¬ 
sition. The most important feature emerging from Fig. 4 
is that the ratio between the relative density and temper¬ 
ature changes associated with the second sound solution, 
have an opposite sign and a much larger value in modulus, 
reflecting that second sound, for temperatures larger than 
the hybridization value is dominated by the fluctuations 
of the density, rather than by the ones of the temperature. 
This important feature is also revealed by the fact that, in 
the same range of temperatures, second sound practically 
exhausts the compressibility sum rule [6] 

(pdP/dplrT 1 = ^ j du>S(q,io)/u , (18) 
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Fig. 3: Isothermal compressibility (pdP/dp^) -1 computed us¬ 
ing the Bogoliubov model (blue curve) and the diagrammatic 
approach of Ref. [12] (green curve). The parameters are chosen 
as in Fig. 1. 


where S(q,u) is the dynamic structure factor, a quantity 
which can be easily calculated within the two-fluid model 
described by Landau’s theory [18]. The fact that the com¬ 
pressibility sum rule (18) is strongly affected by second 
sound is a remarkable feature exhibited by the weakly- 
interacting Bose gas above the hybridization point which 
distinguishes it in a profound way with respect to the case 
of strongly interacting superfluids, like 4 He or the Fermi 
gas at unitarity, where the density fluctuations associated 
with second sound are instead always very small. 

From an experimental point of view this result also re¬ 
veals that second sound is much more easily accessible 
than first sound in dilute Bose gases, being very sensi¬ 
tive to the coupling with a density probe, as confirmed 
by its experimental identification in the experiment of 
Ref. [19]. It is finally worth noticing that the fact that 
second sound in a Bose gas can be easily excited through 
a density probe is not specific to the 3D case. Indeed, 
a similar behavior takes place also in 2D Bose gases [20] 
where its measurement could provide an efficient determi¬ 
nation of the superfluid density, including its discontinuity 
at the Berezinski-Kosterlitz-Thouless transition. Probing 
directly the hybridization mechanism between first and 
second sound is a more difficult experimental task because 
of the extremely low temperatures required to observe 
the coupling between the two modes. The achievement 
of the hybridization could be detected by the observation 
of two sound waves excited by the same density pertur¬ 
bation and propagating with relatively close, but distinct 
velocities. An alternative perspective is provided by the 
study of the hybridization exhibited by the frequencies of 
the discretized modes of a harmonically trapped gas at 
finite temperature, as discussed in [10]. 
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Fig. 4: Ratio ( Sn/n ) x ( t/8t ) for the lower branch (light blue) 
and the upper branch (dark blue). The parameters are chosen 
as in Fig. 1. 
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